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Abstract 


The cylindrical honeycomb sandwich laminate has been analyzied for natural fre- 
quency theoretically and by finite element method. The basic material properties 
Eu, E22, Gi 2, G23, G31, U12, E s , G s are considered to be random variables with known 
statistics. The randomness in the material properties is handled by using a perturba- 
tion technique and the mean and the standared deviation of the natural frequency has 
been obtained for the cylindrical laminate. The boundary conditions are assummed 
to be simply supported at the panel ends. The results are validated by comparision 
with cylindrical composite panel. 

Numerical results for mean and variance of the honeycomb cylinder are also ob- 
tained using FEM with Monte Carlo simulation. However the small size sample is 
not found to yield representative results. 
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Chapter 1 


Introduction 


1.1 Sandwich Construction 

Sandwich with aluminium honeycomb cores or flex cores are used extensively in 
aerospace, land and water transportation and sports industries because of their high 
specific strength and stiffness. Some examples of structures composed of sandwich 
panels are launch vehicles, satellites, helicopter rotor blades, aircraft control surfaces, 
ship hulls, jet engine nacelles and skin, high-speed trains, etc. 

This type of construction consists of two thin facing layers separated by a core 
material. Several types of core shapes and core material have been applied to the 
construction of sandwich structures. The core layer is made of low specific weight 
material like balsa, porous rubber, corrugated metal sheet, metallic and non-metallic 
honeycomb, etc., which may be much less stiff and strong than the face sheets. 

In addition to the possibility of achieving high flexural-stiffness providing a smoother 
aerodynamic surface in a high-speed range, sandwich-type constructions also exhibit 
many properties of exceptional importance for aerospace and civil constructions. 
Among these are: (a) excellent thermal and sound insulation; (b) a longer time 
of exploitation as compared to stiffened-reinforced structures; (c) possibility of being 
designed as to meet very close thermal distortion tolerances such as those required 
for communication satellite antennas and reflectors. 

The most popular core is the honeycomb construction that consists of very thin 
foils in the form of hexagonal cells perpendicular to the facings. Face sheets usually 
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consist of aluminium or fiber-reinforced composite laminates. Face sheets are stiff 
and strong because they carry most of the loads, while lightweight core separates face 
sheets so that a higher bending stiffness of the composite plate can be achieved. 

Regardless of the face sheet and core materials, the sandwich panel is considered 
to be a composite structure because of its inherent in homogeneous and anisotropic 
nature. Sandwich construction is frequently used instead of increasing material thick- 
ness. A sandwich construction provides excellent structural efficiency, for example 
high ratio of strength to weight. Other advantages offered by sandwich construction 
are elimination of welding, superior insulating qualities and design versatility. Even 
though the concept of sandwich construction is not very new, it has only been adopted 
for primary part of structures very recently. This is because there are a variety of 
problem areas to be investigated and overcome when the sandwich construction is 
applied to design of dynamically loaded structures. To enhance the attractiveness 
of sandwich construction, it is thus essential to have a understanding of the local 
strength and dynamic characteristics of individual sandwich panel/beam members. 

Advanced sandwich-type constructions imply the presence of thick orthotropic 
core with bonded anisotropic face sheets that are treated as composite laminates, 
this arrangement presents an opportunity to tailor both the physical and mechani- 
cal properties of the faces by proper selection of laminate materials, their stacking 
sequence and fiber orientation. Suitable selection of fiber orientation and stacking 
sequence can result in substantial improvements of the buckling strength and the 
non-linear response behavior to a variety of load conditions. The transverse shear 
elastic module of the core layer can also be optimized so as to enhance the overall 
response behavior of the sandwich constructions. 

As expected, analytical modelling of sandwich-type panels is much more intricate 
than that of the usual laminated composite structure. In contrast to the case of the 
regular laminated composite structure for which the assumptions are postulated for 
the structure as a whole, in the case of sandwich-type constructions the assumptions 
involving the core layer are different from those associated with the face sheets. More 
over, the analysis of sandwich panels featuring laminated face sheets is much more 
complicated than that with single layered faces. The complexity is due to the presence 
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of three types of asymmetries resulting from the lay-up sequences in the face sheets 
namely: (1) asymmetry with respect to the mid-surface of the face sheets, referred 
to as face asymmetry, inducing face bending-stretching coupling (2) asymmetry with 
respect to the mid-surface of the core, refer to as global asymmetry, which induces 
global bending-stretching coupling, and (3) presence of ply-angles between the prin- 
ciple axis of orthotropy of the face sheets materials and the geometric axes of the 
panel, inducing a structural coupling between stretching and shearing. 

1.2 Literature survey 

1.2.1 Theoretical and Experiment Studies 

Many investigators have carried out noteworthy theoretical and experimental stud- 
ies on linear elastic and nonlinear behavior of aluminium sandwich panels. Paik et 
al. [1] have studied the strength characteristic of aluminium sandwich panels with 
aluminium honeycomb core theoretically and experimentally. A series of strength 
tests have been carried out on aluminium honeycomb-core sandwich panels specimen 
in three point bending, axial compression and lateral crushing. Simplified theories 
have been applied to analyze bending deformation, buckling/ultimate strength and 
crushing strength of honeycomb sandwich panels subjected to the corresponding load 
component. The structural failure characteristics of aluminium sandwich panels have 
been discussed. 

Becker [2] has performed closed-form analysis of the thickness effect of regular hon- 
eycomb core material. In sandwich structure in many cases the effect of a honeycomb 
core in not only to maintain the distance between the face skins, but it also contributes 
to the overall in-plane stiffness. Due to the coupling of the core displacements with 
those of the face sheets the stiffness contribution is not simply proportional to its 
total thickness but is a nonlinear function of the core thickness which is the so-called 
core thickness effect. The closed-from analysis is given for the effective in-plane core 
stiffnesses including the thickness effect. Based on an approximate representation for 
the displacement field within the core cell walls, the effective core stiffness have been 
derived by energy considerations. The results have been validated by comparing with 
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finite element analysis [FEM]. 

Fatt and Park. [3] have studied experimentally a three stage, progressive damage 
model for the perforation of a honeycomb sandwich plate subjected to normal im- 
pact by blunt and spherical projectiles. The residual velocities in three stages have 
been found from energy balances between each stage. The plastic work dissipated 
in deformation and fracture at each stage have been approximated from the solution 
of the previous stage. Shear forces have been transmitted in the bond between face 
sheets and honeycomb core. The predicted ballistic limits have been found to be 
within 5% of the measured ballistic limits for sandwich plates perforated by the blunt 
and spherical projectiles. The lateral extent and deformation of the damaged face 
sheets have been also found to be with 34 and 51% of limited test results. They have 
found that most of the work is dissipated in perforation of the bottom face sheet dur- 
ing impact. Li [4] has obtained stress concentration and local behavior of prestressed 
composite laminates/sandwich shells with overhangs in cylindrical bending. The lam- 
inates have been considered as cross-ply /angle-ply multilayer shells with overhangs in 
cylindrical bending caused by transversely concentrated line-load and in-plane uni- 
form load. The high stress gradient, wich could not be measured by experimentally, 
Li has determined theoretically around the load area and sharply support regions. 
The experimental results have been compared with FEM. Li has used FEM package 
operating system MAEC. 

Glass et al. [5] have analysed Graphite/epoxy honeycomb core sandwich perme- 
ability under mechanical loads. The core material for the test specimens were either 
Hexcel HRP-3/16-8.0 or DuPont Korex-1/8-4.5 and had nominally 1.27 cm thick- 
ness. The facesheets were made of Hercules’ AS4/8552 graphite/epoxy composite 
with nominally 0.15 cm thickness. The sandwich specimens were studied during both 
static and dynamic shear loads. The permeability was measured as the rate of air 
flow through the core from a 2.54 cm diameter circular area of the core exposed to 
an air pressure of 68.9 kPa. During both static and dynamic testing, the Korex core 
experienced sudden increase in core permeability corresponding to the core catas- 
trophic failure while the HRP core generally experienced a gradual increase in the 
permeability prior to a bond line failure. 
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Petras and Sutcliffe. [6] have presented failure mode maps for honeycomb sand- 
wich panels. Theoretical models using honeycomb mechanics and classical beam 
theory are described. A failure mode map for loading under 3-point bending have 
been constructed. The dependence of failure mode and load on the ratio of skin 
thickness to span length and honeycomb relative density has been obtained. Beam 
specimens are tested in 3-point bending. The experimental data agree satisfactorily 
with the theoretical predication. Muc and Zuchara. [7] have performed the analysis 
of a thin- walled sandwich laminate composite faces subject to axial compression. The 
analytical results have been compared with FEM NISA II package. 

Cho and Averill. [8] have studied the sandwich construction with three-dimensional 
finite element based on first-order zig-zag sublaminate approximations. The in-plane 
displacement fields in each sublaminate are assumed to be piecewise linear function 
and vary in a zig-zag fashion through the thickness of the sublaminate, the zig-zag 
function have been evaluated by enforcing the continuity of transverse shear stresses 
at layer interfaces. The in-plane displacement fields have been assumed to vary lin- 
early through the thickness. The transverse normal strain predication have improved 
by assuming a constant variation of transverse normal stress in each sublaminate. 
In the computational model, each finite element represents one sublaminate. The 
finite element has been developed with the topology of an eight-noded brick, allowing 
the thickness of the shell to be discretized into several elements, or sublaminates, 
where each sublaminate can contain more than one physical layer. Each node has 
five engineering degrees of freedom, three translation and two rotation. Shafizadeh 
and Seferis. [9] have performed scaling of honeycomb compressive yield stresses. A 
honeycomb-scaling factor and geometric end constraint have been found to relate the 
rings and core through the relative yield stresses and the physical dimensions and 
properties of the honeycomb materials. The compressive properties of the Nomex 
rings have been also investigated using both a model and a commercial phenolic hon- 
eycomb dip resin. The Nomex rings manufactured from the model resin were found 
to have higher compressive properties. These properties were attributed to the higher 
fracture toughness of the resin, and both resins were found to accurately scale from 
rings to core with the honeycomb-scaling factor. 
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Guo and Gibson. [10] have analysed the behaviour of intact and damaged honey- 
combs using finite elements. The Young’s moduli, the elastic buckling strength and 
the plastic collapse strength of regular honeycomb with defects consisting of miss- 
ing cells in the structure have been analyzed using the finite element method. The 
behavior of intact honeycombs have been found to be consistent with those of pre- 
vious analyses. The effect of single isolated defects of varying sizes and the effect 
of the separation distance between two defects on the elastic and plastic behaviors 
has been studied. Single isolated defects reduce the modulus and strength. The 
elastic strength of honeycomb with a defect normalized by the intact strength de- 
crease directly with the ratio of the minimum net cross-sectional area normalized 
by the intact cross-sectional area. The plastic collapse strength of honeycomb with 
a defect normalized by the intact strength decreases less rapidly than the ratio of 
the minimum net cross-sectional area normalized by the intact cross-sectional area. 
Burton and Noor. [11] have obtained the assessment of continuum models for sand- 
wich panel honeycomb cores. Finite element models have been used for predicting 
the free-vibration response of infinitely long and rectangular sandwich panels. The 
panels considered have square-cell honeycomb core and simply supported edges. The 
sandwich core and face sheets have been modelled as three-dimensional solid elements 
and two-dimensional plate elements. The predications of the finite element models 
have been compared with those obtained by using higher-order sandwich theory for 
panels with the core replaced by an effective continuum. 

Palazotto et al. [12] has analysed the response of composite sanwich plate to 
low-velocity impact predicted by a displacement-based, plate bending, finite element 
algorithm. Fifth order Hermitian interpolation allows three-dimensional equilibrium 
integration for transverse stress calculation to be carried out symbolically on the 
interpolation function so that transverse stress within the elements are expressed di- 
rectly in terms of nodal quantities. Nomex honeycomb sandwich cores have been 
modeled using an elastic-plastic foundation and contact loading have been simulated 
by Hertzian pressure distribution for which the contact radius have been determined 
iteratively. Damage predication by failure criteria and damage progression via stiff- 
ness reduction have been employed. Vaidya et al. [13] have studied the impact 
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response of integrated hollow core sandwich composite panels. In comparison with 
traditional foam and honeycomb cores, the integrated space core provides a means to 
route wires/rods, embed electronic assemblies, and store fuel and fireretardant foam, 
among other conceivable benefits. The low velocity impact response of innovative inte- 
grated sandwich core composite have been investigated. Three thickness of integrated 
and functionality-embedded E-glass/epoxy sandwich core have been considered. The 
low-velocity impact results indicated that the hallow and functionality-embedded in- 
tegrated core suffered a localized damage state limited to a system of core members in 
the vicinity of the impact. Results have been compared with experimental findings. 
The specimens were loaded in compression at a loading rate of 1.27 mm/min and 
taken to failure. The low velocity-impact test have been conducted in a Dynatup 
8120 impact-testing machine equipped with a load cell capacity of 7700 kg. The sam- 
ples were clamped using pneumatic assit, such that a circular opening of 76.2 mm 
diameter have exposed to the impactor. The impactor had 12.5 mm hemispherical 
tup. The force- time, energy-time and load-deflection response of the samples have 
been measured using the instrumented feature of the Dynatup machine. 

Meraghni et al. [14] have studied the mechanical behaviour of sandwich panel 
analytically performed a finite element analysis , an and experimentally. The ef- 
ficiency of a new analytical method for sandwich panels core has been presented. 
The separation of the skins by the core increase the inertia of the sandwich panel, 
the flexure and shear stiffness. Improved mechanical behaviour has been obtained 
with a little increase in weight, producing an efficient structure to resist bending and 
buckling loads. The approaches (theoretical and experimental) have been used to 
determine elastic properties and ultimate stress. Daniel and Abot [15] have studied 
the flexural behaviour of composite sandwich beams and compared the results with 
predications of theoretical models. Sandwich beams have been fabricated by bond- 
ing unidirectional carbon/epoxy laminate face sheets to aluminium honeycomb cores 
with an adhesive film. All constituent materials (composite laminates, adhesive and 
core) have been characterized independently. Special techniques have been developed 
to pervent premature failure under the loading pins and to ensure failure in the test 
section. Sandwich beams have been tested under four-point and three-point bending. 
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Strain to failure in the face sheets have been recorded with strain gauges, and beam 
deflection and strain in the honeycomb core have been recorded by using moire tech- 
niques. The beam face sheets exhibited a softening non-linearity on the compression 
side and a stiffening non-linearity on the tension side. 

Thomsen [16] has presented high-order theory for the analysis of multi-layer sand- 
wich panels. Multi-layer composte assemblies are characterized by having interchang- 
ing layers of high density and high stiffness separated by complaint interface layers. 
A characteristic feature of these assemblies is that severe and highly localized stress 
concentrations appear under certain conditions. The quantitative understanding of 
these phenomena involves a detailed description of the complicated mechanical inter- 
action of the stiff layers through the compliant interface layers. A high-order plate 
theory formulation, which includes the through-thickness flexibilty, and which inher- 
ently includes a description of the global response as well as the local responses of 
the individual layers of an arbitrary, multi-layer structural plate assembly. Wang 
and Yang [17] have performed experimental investigation of the damping behaviour 
of laminated honeycomb cantilever beams with fine solder balls enclosed in the cells 
as dampers. Attenuation have been achieved through the exchange of momentum 
between the balls and the beam. The effect of the mass ratio ( the number of balls 
to be enclosed in each honeycomb cell) have been studied by varying the number of 
balls used. The damping variation have been found to be effective in reducing the 
amplitude without significantly shifting the natural frequency of the cantilever. 

1.2.2 Imperfection in System Parameters 

As is well known, many uncertainties exist during the process of fabrication measure- 
ment and manufacture of structures. A composite laminate is a stack of layers of 
fiber-reinforced laminae. The fiber-reinforced laminae are made of fibers and matrix 
that are of two different materials. The way in which the fiber and matrix materials 
are assembled to make a lamina, as well as lay-up and curing of laminae, are compli- 
cated processes and may involve a lot of uncertain factors. Therefore, the material 
properties and geometry of composite laminate are random in nature. 

In many engineering applications, the structures may be subjected to different ran- 
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dom loads and the corresponding failure modes may be different as well. If the loads 
are not applied simultaneously, there may be some stochastic correlation between the 
different loadings. This leads to the necessity to take account of the uncertainties 
of loads and some other design parameters in the reliability analysis. More effort is 
need in developing analysis techniques to obtain the combined effect of uncertainties 
in applied loads, material properties, fiber orientations and plate thickness on the 
reliability of composite laminates. This attaines special significance if highly reliable 
composite structures are to be designed. 

It’s well known that discrepancies between analytical and experimental predictions 
of response of sandwich panels and shells are mainly due to the presence of unavoidable 
geometric and material imperfections. For this reason, in recent years the study of the 
implications of the imperfection on the load carrying capacity of sandwich structures 
and its influence upon the structural design has been very important. 

Chen et al. [18] have performed a finite element study of the effect of holes and 
rigid inclusion on the elastic modulus and yield strength of regular honeycombs under 
biaxial loading. The focus honeycomb were weakened by a small degree of geometrical 
imperfection, such as a random distribution of fractured cell walls. It has been found 
that the strength of an imperfect honeycomb is relatively insensitive to the presence of 
holes and inclusions, consistent with recent experimental observations on commercial 
aluminium alloy foams. 

A problem of considerable practical importance towards accurate analysis and en- 
hancement of load carrying capacity of sandwich curved panels consists of a better 
understanding of the possible implication of material imperfections, and implemen- 
tation of the yielding mechanism and the attenuation of imperfection sensitivity to 
unavoidable geometric and material imperfections. However, in spite of the evident 
importance of the problem, to the author’s knowledge, no results related with these 
issues have been reported in the specialized literature. 
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1.3 Present Work 


Uncertainity in structural dynamics mainly arise from the uncertainly in the struc- 
tural properties. The problem is statistical in nature and is due, for example, to 
the stiffness or damping fluctuation caused by random variations in material prop- 
erties, randomness in boundary conditions, and variations caused by manufacturing 
and assembly techniques. 

The aim of the present study is to investigate the dynamic characteristics of the 
honeycomb-cored sandwich cylinders under material imperfection. The material prop- 
erties are modeled as random in nature. The system equations involve partial differ- 
ential equations whose coefficients are random function of space. 

This modelling of the material properties as random variables results in random 
eigenvalues, eigenvectors, and response of system in question. In general, it is possible 
to obtain close from expression for the response statistical characteristics for only 
specific simple cases of geometry, boundary conditions and loadings. Cases for which 
analytical solutions are not available, some numerical approach has to be used. The 
numerical methods include the Monte Carlo simulation and finite element method. 
These, through computational intensive and time consuming, can be applied to any 
problem in this area. The present study has used the two in combination to generate 
the second order statistics of the response of the sandwich cylinder. The simply 
supported edge conditions have been considered for the cylinder. The effect of shear 
deformation has been incorporated in the study. 
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Chapter 2 


Theory 


The behaviour of the sandwich panel may be studied by considering the behaviour 
of the honeycomb core and the facing lamina separately. The performance of the 
complete sandwich is obtained by combining the two portion together. 

2.1 The In- Plane Properties of Honeycomb 

The honeycomb portion is not a homogeneous section. The behaviour of the thin 
walled cells under loading is not similar to a solid lamina. The thin walls under load 
acts as slender beams. The performance of the honeycomb is being examined under 
inplane loading of uniaxial and biaxial types. 

2.1.1 Uniaxial Loading 

The in-plane stiffness and strength, that is, (Figure 2.1) those in the X\ - X 2 plane 
of hexagonal honeycomb panel are the lowest because stresses in this plane make the 
cell walls bend. The out-of-plane stiffness and strength those in X 3 direction are 
much larger because they require the axial extension or compression of the cell walls. 
The out-of-plane analysis is carried out separately as it gives the additional stiffnesses 
which are needed for the design of honeycomb cores in sandwich panels and for the 
complete description of the behaviour of natural honeycomb. Figure 2.1 shows a unit 
cell of the hexagonal core. The figure shows the undeformed cell with dimensions, h 
is the cell wall length in X 2 direction, l is the length of the inclined wall, t is the wall 
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Figure 2.1: Unit cell of hexagonal honeycomb 
(a) undeformed honeycomb, (b) and (c) bending due to loads in the X x and X 2 directions 

thickness and 6 gives the angle of inclination. It is assumed that the general cell may 
have different thickness for different walls. We wish to analyse its response to load 
applied in the X x - X 2 plane, if the hexagonal is regular, that is, the side are equal 
and the angle are all 120° and the cell walls are of the same thickness, then the in- 
plane properties are isotropic: they do not depend on direction. Such a structure has 
two independent equivalent elastic moduli a Young’s modulus E and a shear modulus 
G, for instance. But when the hexagon is irregular or the cell walls in one direction 
are thicker than those in the others, the properties are anisotropic and a complete 
description of the in-plane properties requires five modulii, say E x , E 2 , G X2 , u u and v 2x . 
When a honeycomb, loaded in the X x or the X 2 direction, deforms in linear-elastic 
way, the cell walls bend. The response is conveniently described by the above five 
moduli, but the five are not independent. The reciprocal relation reduces the number 
of independent moduli to four. They are the same in tension as in compression. 


12 




Consider a unit cell, as shown in figure 2.1a, to include two inclined walls and 
one wall parallel to X 2 axis. A force component P giving rise to stress <J\ parallel 
to Xi causes one set of cell walls-those of length l - to bend as shown in Figure(2.1 
b). Equilibrium requires that the component of forces C parallel to X 2 be zero. The 
moment M acting to bend the cell wall, which we treat as a beam of length l, thickness 
t, depth b, and Young’s modulus E s , is from moment equilibrium 
Plsind 


M = 


where, from force equilibrium in X\ direction with b as the height of the cell 


P = <7i (h + Isind ) b 


( 2 . 1 ) 


( 2 . 2 ) 


From standard beam theory, modelling the wall as a cantilever beam, the wall deflects 
by: 


5 = 


Pl z sin9 

12EJ 


(2.3) 


where I is the second moment of inertia of the cell wall (J = bP / 12 for a wall of 
uniform thickness t ). A component of the deflection SsinO is parallel to the Xi axis, 
giving a strain: 

5 sinQ , , 

€l = 7 —— 7T ( 2 -4) 


Using eq (2.3) 


Substitute eq (2.2) in eq (2.5) 


lcos9 

Pl 2 sin 2 9 
61 ~ 12 E s Icos9 

( h + IsinO) bPsirPd 


(2.5) 


( 2 . 6 ) 


12E s Icos9 

The Young’s modulus parallel to X\ is just £j = cq/ei, giving the equivalent modulus 

_ < 7112^/0050 

1 <7i (h + lsin9)bl 2 sin 2 9 

For a cell with uniform wall thickness t 

3 E,cos9 


MDV 


(2.7) 


(j + sin9 ) sin 2 9 

Proceeding as above for the X 2 direction, the force acting on the cell wall of length, 
l and depth, b is given by 


W = a 2 lbcos9 
13 



When <72 is the stress parallel to X 2 . The bending moment acting on the wall is 


M 


Wlcosd 


The wall taken as cantilever beam, deflects by: 

Wl 3 cos0 
12 E S I 

Of this, a component 8cos9 is parallel to X 2 axis, giving a strain: 

ScosQ 

h + IsinQ 

Substituting for S from eq (2.9) and for W in-terms of S 2 gives 

o 2 lbcosQl 3 cos 3 0 


62 


After simplify 


€2 = 


12 E s I(h + IsinQ ) 
o 2 bl 4 cos 3 Q 


12 E s I(h + IsinQ ) ' 

The Young’s modulus parallel to X 2 is simply E 2 = cr 2 /e 2 , giving: 

„ _ <j 2 12 j E' s J(/i + IsinQ) 

2 a 2 bl 4 cos 3 Q 

For a beam of constant thickness, thesis becomes 

't\ 3 E s (j + sind'j 


Eo 


G)‘ 


cos 3 6 


( 2 . 8 ) 


(2.9) 


( 2 . 10 ) 


( 2 . 11 ) 


( 2 . 12 ) 


The poisson’s ratios are calculated by taking the negative ratio of the strains 
normal to, and parallel to, the loading direction. We find for loading in Xi direction: 


^12 = 

and for loading in X 2 direction: 

^21 = 


- 

cos 2 6 

€l 

(j + sind\ sinQ 


£1 _ 

( j + sinOj sin6 

62 

cos 2 6 


(2.13) 


(2.14) 


The shear modulus, G 12 can be calculated using Figure 2.2. Because of symmetry 
there is no relative motion of the points A, B and C when the honeycomb is sheared; 
the shear deflection u s is entirely due to the bending of beam BD and its rotation 

14 




A c 

fb - ) 


Figure 2.2: Cell deformation by cell wall bending and rotation due to shear 

(a) the undeformed honeycomb (b) loads, moments, displacements and rotations due to shear. 
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(through the angle 4>) about the point B. The force are shown in the figure. Summing 
moments at B gives the bending moment applied to the members AB and BC: 


M = 


Fh 

4 


Then, using the standard result for cantilever beam 


6 = Ml 2 /6E S I 


(2.15) 


The angle of rotation is: 


4>' 


Fhl 
24 EJ 


The shearing deflection u s of the point D with respect to B is: 


1„ F fh ' 

“* = 2 #+ 3^/U 


Fh 2 
48 EJ 


(i l + 2 h) 


The shear strain 7 is given by: 


2u s 


Fh 2 (l + 2 h) 


^ (h + IsinO) 2AE S I ( h + IsinO ) ^ ^ 

The shear stress is r = F/2lbcos9. The shear modulus can be obtained as G 12 — T lh 


G 


12 


12E a I ( h + IsinO ) 
MFc^sO (1 + 2h) 


For a cell with uniform wall thickness 

G12 = 


t \ 3 E s (j + sinO) 

O (s ) 2 (1 + f) cose 


(2-17) 


2.1.2 Biaxial Loading 

The biaxial in-plane loading of the honeycomb can be visualised with the help of 
Figure 2 . 1 . The principal stresses are assumed to be oi and <r 2 and P and W are 
the forces in X x and X 2 directions. The low in-plane stiffness of honeycomb is due to 
cell-wall bending. In uniaxial loading the bending displacements are so much larger 
than the axial stretching or compression of the walls that these can be neglected. For 
the biaxial loading the effects of the axial strains are incorporated below. 

Consider the strains due to bending first. Referring to Figure 2 . 1 , the force com- 
ponents P in the X x direction causes the cell walls of length l to bend in one direction. 
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In the biaxial loading a force W in the X 2 direction in superimposed. The net dis- 
placement is calculated by the method of section 2.1.1 and from this the strain e x and 
e 2 are obtained. 

£l= (a . 18) 
€2 = ^{<*2 - * 2 lOl ) 

The stresses cq and <r 2 introduce an axial force F x on the walls of length l, and an 
axial force F 2 on the walls of length h. 


Fi = PcosQ + WsinQ F 2 = 2 W 


The additional strain in X x direction is: 

PcosQ + WsinQ 
Cl = bt. E s 


(2.18a) 


In the X 2 direction the additional strain is: 

1 f ( PcosQ + W sinQ)lsinQ + 2Wh 
62 btE s \ (h + IsmQ ) 

substituting for P and W from the equations in section 2.1.1 gives the additional 

strains in terms of t, l, h and Q. Adding these to the contribution from bending gives 

the full elastic consecutive equation for in-plane deformation: 


(2.186) 


e i = E ^( a i ~ V 1202 ) + £^ry {cq (j + sinQ^j cosQ + <j 2 sm0cos^| 
C 2 = ~ + {^l sinQcosQ + ^2 } 


(2.19) 


2.2 Out-of-plane properties of Honeycombs 


Here we note that the function of the honeycomb core is to carry normal and shear 
loads in planes containing the axis of the hexagonal prisms in the Xz direction. When 
loaded in this direction the cell walls are extended or compressed rather than bent 
and the moduli can be expressed as Ref[22] 


F _ f - f ML + 2 \ 

3 s l \2(h/l + smQ)cos9 J 


( 2 . 20 ) 


The strain energy can be calculated using Figure (2.3) . A set of displacements 
are assumed which satisfy the external boundary conditions. Consider the uniform 
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Figure 2.3: Honeycomb loading face to the normal 


shear strain 7 13 and associated shear stress r 13 acting on the normal face to X 3 in 
the Xi direction. The elastic strain energy is stored in the shear displacements in the 
cell walls. The bending stiffnesses and the energies associated with bending are much 
smaller. The shear strains in wall a, b and c are 


7a = 0 

7 6 = 7i3 cosd 
7c = 7i3 cosd 


A lower bound for moduli is Ref[22] 


n , GjCosQ (t\ 
13 — (h/l + sind ) \l) 


( 2 . 21 ) 


The stress distributions satisfy equilibrium at each point and are in equilibrium 
with the external loads. Consider first loading in the X x direction. An external stress 
r 13 induces the set of shear stresses r a , n and r c in the three cell walls. By symmetry 
and as the wall a is loaded in simple bending it does not carry significant load are 


T a — 0 


And 


n = T c 

2ri 3 (h + lsinQ)cos6 = 2Tbtlcos9 


(2.21a) 
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The equilibrium equation, gives a lower bound 

G s cos9 


G 13 > 


(h/l + sin9 ) 


(D 


(2.216) 


For regular hexagons it reduces to 

G13 = 0.577G, Q) (2.21c) 

Due to an external shear strain 7 23 in the X 2 direction the strains in the cell walls 


axe 


7 a = 723 

7& = 723 sinO 
7c = 723 sinO 


And a lower bound the moduli is 


1 G s (h/l + 2 sin 2 9) ft 
G23 < pr 


G) 


2 (h/l + sin9)cos9 \lj (2.22) 

Consider loading in the X 2 direction an external stress r 2 3 induces the set of shear 

stresses r a , r b and r c in the three cell walls. 

U = t c 

Equilibrium in the X 3 direction 

r a = n + r c = 2 r b 

and equilibrium with the external stress gives 


2r 23 l(h + lsin9)cos9 = 2 r b tlsin9 + r a th 

l 


G 23 > 


n = T 23 COs9j 
G s (h/l + sin9) ft 


G) 


(2.22a) 


(2.226) 


(2/i/J + l)cos# 

In this case the two bound do not coincide for a general, anisotropic honeycomb, even 
though they are close. But for regular hexagons it reduces to 

G 23 = 0.577G* Q) (2.22c) 

The two Poisson’s ratios 1/31 and u 32 are simply equal to those for the solid 

^31 = ^32 = (2.22 d) 

The Possion’s ratio’s v X3 and v 23 are found from the reciprocal relation 

0. 


Ei ~ 
u l3 = —v s « 0 , 
■E3 


E 2 

V 23 = — v 8 


(2.22e) 
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2.3 Stress-Strain relation for honeycomb 

The three-dimensional constitutive relations for stress and strain for the honeycomb 
panel can be expressed as 


where 


en E ^ <Jl1 z/l2<7 22) + +<22202 — — 3033 

1 1 El 
622 = Y (~ U2lUl1 + ^ 22 ) + cr n a 2 + a 22 a 3 - 


_ i , 

633 ~ E — l / 32 (j 2 2 + (7 33 ) 


712 

J_ 

723 

_ 1 _ 

731 


G 


12 


•7-12 

<^23 

723 

^31 

731 


a 1 
a 2 


(h/l+sinQ)cos9 

sinOcosQ 

W T) 


E 3 (t/l)(h/l+sin9) 


In matrix notation eq (2.23) to (2.28) can be expressed as 


Eo 


(2.23) 

(2.24) 

(2.25) 

(2.26) 

(2.27) 

(2.28) 


e ll 


ST + a i 

-S+“2 

Ez 

0 

0 

0 


' N 

o-n 

e 22 


“If + “2 

i + 0 3 

__ ^32 

Ez 

0 

0 

0 


<722 

e 33 

> _ 

— iOi 

Ex 

__ ^21 

E 2 

X 

Ez 

0 

0 

0 


<733 

723 


0 

0 

0 

1 

G 23 

0 

0 

< 

► 

723 

731 


0 

0 

0 

0 

1 

G 31 

0 


731 

712 , 


0 

0 

0 

0 

0 

1 

G 12 - 


s r 12 , 


(4 = [5] {o} 

Inverting the above relation 
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(2.30) 


\ 

Cll 


k 

P 

k 

p 

k 

/3 

0 

0 

0 


f > 

Cll 

<?22 


k 

(3 

k 

p 

fie. 

/3 

0 

0 

0 


622 

033 

y _ 

k 

P 

ks. 

p 

k 

P 

0 

0 

0 

< 

£33 

► 

723 


0 

0 

0 

G 23 

0 

0 


723 

T 31 


0 

0 

0 

0 

G 31 

1 

0 


732 

, T 12 , 


_ 0 

0 

0 

0 

0 

G l2 _ 


. 712 , 


where for the sake of compactness j3 denotes |/ 3 ] 


m 


bi- 

h 
b 3 
b 4 

h 

b 6 

b 7 
&8 
&9 


1 I a 3 | a i | a 3 Qi ^ 32^23 01 ^ 32^23 

S 1 S 2 S 3 ' Si S 3 " r S 2 S 3 ‘ S 3 S 1 S 2 S 3 E 2 E 3 

_f 2 _ 4 . M 12 4 . £2J31 _ J^lEVL. _!_ Q2 g iag32. _ EtiElZEZZ 
E 3 E 1 E 3 E 2 E 3 E 1 E 2 E 3 E 1 E 3 E\ E 2 E 3 

-1- ^2^31^12 ___ t/12^23^31. __ .*43,£31 £3J£JL3.Q3. 

' E 1 E 2 E 3 E 1 E 2 E 3 E 1 E 2 E 3 E 1 E 3 

1 _ £32S23. 

E 2 E 3 E 3 S 2 E 3 

_£ 21 _ 4 _ £23 £31 

S 2 S 3 S 2 S 3 

. ryoi 1 tvoi 


_ QTL. _L_ 

S 3 ^ 


93E1 1 4 . ¥21E21 4 . _£2JL 
S 3 S 2 S 3 S 3 S 


^31 1 Q 3^31 

S 3 S 2 S 3 


ju /3 

_ ^ 2 . j_ _J42_ 4 . mS32. 

S 3 Si S 3 Si S 3 

1 | CL\ ^31 ^13 

E 1 E 3 E 3 E 1 E 3 

E 2 . 2 4. ^32 °1 _ ^ 31 Q 2 . 1 ‘ • j.im 

Ei S 3 S 3 S 3 Si S 3 

_ q.2S23. 4 . £ 12 S 22 4 . -SI3— 4 . aa£l3 
S 2 ^ S 1 S 2 ^ S 1 S 2 ^ Si 

^ 1 3Q J2 4 . £2L£l3. 
Ei Ei S 2 


^23 | ^23 Ql 

Si s 2 s 2 

_JL 

Si 

These 


M = [ 0 ] (4 


£ 21 * 42 . 

S 1 S 2 


where 5' 


<3 


[Qij] — 


k 

p 

k 

P 

k 

/3 

0 

0 

0 

k 

J 9 

k 

P 

k 

P 

0 

0 

0 

k 

P 

k 

(3 

k 

/S 

0 

0 

0 

0 

0 

0 

G 23 

0 

0 

0 

0 

0 

0 

G 31 

0 

0 

0 

0 

0 

0 

G 12 


or 


(2.31) 
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where 


A.i A-2 -43 0 0 0 

A 4 A.5 A g 0 0 0 

Aj As A9 0 0 0 


0 0 0 Go 


0 0 


0 0 0 0 G zl 0 

0 0 0 0 0 G\ 


A\ — 


Ao — 


A 3 = 


a ~ h. a - h. a -h. a - h. 

6 \py 7 ifl’ 8 w 9 m 

The laminate theory is valid for composite materials which are laminated in thin 
layers. Core materials constituting the sandwich panel have a high thickness. Hence 
sandwich panel can’t be studied by using the laminate theory. In this theory it is 
possible to determine the equivalent rigidities. Laminate theory can be modified 
to take into account the transverse shear. Laminate approach also can be used for 
sandwich panel. A regular unit cell is consider with cell parameters l and h. X 3 axis 
is directed vertically upwards. The honeycomb walls is considered to be one solid 
layer of laminate. The wall thickness in the three inclined directions of the hexagonal 
shape are used to evaluate the equivalent geometric and material parameters of this 
lamina. The transformation matrix around A 3 -axis is defined by: 


cos 2 6 sin 2 9 0 0 

sin 2 6 cos 2 6 0 0 


2cos9sin9 
—2 cos9sin9 


0 0 0 cos9 —sin9 0 

0 0 0 sin9 cos9 0 

—cos9sin9 cos9sin9 0 0 0 cos 2 9 — sin 2 9 


(2.32) 
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[T]- 1 


cos 2 9 
sin 2 9 
0 


sin 2 9 0 0 0 

cos 2 9 0 0 0 

0 10 0 


0 0 0 cos9 sin9 

0 0 0 —sin6 cos9 

cos9sin9 —cos9sin9 0 0 0 

Where [R] is the Router matrix defined by: 


—2cos9sin9 

2cos9sin9 

0 

0 

0 

cos 2 9 — sin 2 9 


(2.33) 


and 


[R} = 


1 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

2 

0 

0 

0 

0 

0 

0 

2 

0 

0 

0 

0 

0 

0 

2 



1 0 0 0 0 0 

0 1 0 0 0 0 

0 0 1 0 0 0 

0 0 0 1/2 0 0 

0 0 0 0 1/2 0 

0 0 0 0 0 1/2 


The relationship between [Qij] and [Q lJ ] is: 


[Qij] = where ij=l to 6 



Qi i 

Q 12 

Ql3 

Q 14 

Q 15 

Ql6 

Q 21 

Q 22 

Q23 

Q24 

Q25 

<^26 

Q31 

Q 32 

Q 33 

Q 34 

Q35 

<?36 

Q41 

Q 42 

Q43 

Q44 

Q45 

^46 

Q51 

Q 52 

Q53 

Q 54 

Q 55 

<3 56 

^61 

Q 62 

Q 63 

^64 

^65 

<366 


(2.34) 


(2.35) 
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where the elements are expressed as 

Qu = A l cos 4 d + A 4 cos 2 9sin 2 9 + A 2 cos 2 9sin 2 9 + A 5 sin 4 0 + AG X2 cos 2 9sin 2 9 
Qn = A x sin 2 9cos 2 9 + A 4 sm 4 6 + A 2 cos 4 9 + A b sin 2 9cos 2 9 - 4 G l2 cos 2 0sin 2 0 
Q13 = A$cos 2 9 + A b sin 2 9 
Ql4 = 0 
Ql5 — 0 

Qie = A l cos 3 9sin9 + A 4 cos9sin 3 0 - A 2 cos 3 9sin9 - A b cos9sin 3 9 
-2 G l2 cos 3 9sin9 + 2Gi 2 cos9sin 3 9 

Q21 = A lC os 2 9sin 2 9 + A 4 cos 4 9 + A 2 sin 4 9 + A 5 cos 2 0sin 2 9 - 4 G 12 cos 2 9sin 2 0 
Q22 = Aisin 4 9 + A 4 cos 2 9sin 2 9 + A 2 sin 2 9cos 2 9 + A b cos 4 9 + 4 G l2 cos 2 9sin 2 6 
Q23 = A z sin 2 9 + A 6 cos 2 9 
Q 24 = 0 
Q25 = 0 

Q 26 = A 1 cos9sin 3 9 + A 4 cos 3 9sin9 - A 2 sin 3 9cos9 - A 5 cos 3 9sin6 
+2 G 12 cos 3 0sin0 — 2Gi 2 cos9sin 3 9 
Q31 = Ajcos 2 9 + A$sin 2 9 
Q32 = A r sin 2 9 + A 8 cos 2 9 
Q33 — A$ 

Q34 — 0 
Q35 = 0 

Q 36 — A 7 cos9sin9 — A 8 cos9sin9 
Q 4x — 0 
Q42 = 0 
Q43 = 0 

Q44 = G 23 cos 2 9 + Gzisin 2 9 

Q45 — —G 28 cos9sin9 + G 8 icos9sin9 

Q46 = 0 

Qsi = 0 

Q52 — 0 

Q33 = 0 

Q54 = —G 28 cos9sin9 + G 8x cos9sin6 
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$55 = G 23 sin 2 8 + G 3 icos 2 8 

$56 = 0 

$61 Aicos OsinO A 4 cos 3 9 sinQ + A2cos9sin 3 0 — A3 cosOsin 3 6 
— 2 Gi 2 cos 3 9 sin 9 + 2Gi2Cos8sin 3 9 

$62 = A x cos0sin 3 6 - A 4 cos6sin 3 0 + A 2 cos 3 9 sin 9 - A 5 cos 3 6sin8 
+2 G X2 cos 3 9 sin 9 — 2G X 2sin 3 6cos9 
$63 = AzcosOsinO — A§cos 9 sin 8 

$64 = 0 

$65 = 0 

$66 = Aicos 2 6 sin 2 9 — A 4 cos 2 6sin 2 6 — A2cos 2 0sin 2 9 + A5Cos 2 9sin 2 9 
+G X2 cos 4 9 — 2Gi2Cos 2 9sin 2 9 + G X 2sin 4 9 
The stress-strain relation may be put as 


( > 
O-ll 


$11 

$12 

$13 

0 

0 

$16 


' \ 

£n 

022 


$21 

$22 

$23 

0 

0 

$26 


e 2 2 

033 


$31 

$32 

$33 

0 

0 

$36 

< 

£33 

723 


0 

0 

0 

$44 

$45 

0 


723 

731 


0 

0 

0 

$54 

$55 

0 


731 

T X2 

V > 


. $61 

$62 

$63 

0 

0 

$66 _ 


. 712 , 


Substituting eq. (2.22d) and eq. (2.22e) in the above relation 



$ 

13 = 

$23 = 

$63 

= 0. 



/ ■N 

0 11 


$11 

$12 

0 

0 

0 

$16 


' \ 

£ll 

022 


$21 

$22 

0 

0 

0 

$26 


£22 

033 

^ __ 

$31 

$32 

$33 

0 

0 

$36 

< 

£33 

\ 

r 23 


0 

0 

0 

$44 

$45 

0 


/ 

723 

731 


0 

0 

0 

$54 

$55 

0 


731 

. r 12 , 


. $61 

$62 

0 

0 

0 

$66 


. 712 , 


(2.36) 


2.4 Shear Deformable Curved Sandwich Panel 


Figure 2.4 show the orthogonal curvilinear coordinates or shell coordinates for the 
sandwich curved panel. The axis carry dual naming a£Xi,X 2 ,X 3 and X,Y,Z. These 
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axe used interchangeably in the analysis. The X and Y axis are along the lines of 
curvature on the mid-surface z=0, and Z-curves are straight lines perpendicular to 
the surface z=0. For cylindrical and spherical shells the lines of principal curvature 
coincide with the coordinate lines. The values of the principal radii of curvature of 
the middle surface are denoted by Ry and R 2 . 

The position vector of a point on the middle surface is denoted by r, and the 
position of a point at distance z from the middle surface is denoted by R. The distance 
ds between points (x,y,0) and (x+dx,y+dy,0) is determined by [19,20] 


{dsf 


= dr. dr 

= cti 2 (dx ) 2 + a 2 2 {dy ) 2 


(2.37) 


where dr = tydx + t 2 dy , the vector ty and t 2 (t t = |£), i=l,2 are tangent to the Xy 
and x 2 coordinate lines, and ay and a 2 are the surface matrices 


a\ = Mi, 


& 2 ~ ^ 2-^2 


(2.38) 


The distance dS between points(x,y,z) and (x+dx, y+dy, z+dz) is giving by 


( dS ) 2 = dR.dR 

= Ll(dx) + Lj(dy) + L 2 (dz) 


(2.39) 


where dR = ^dx + ^ jdy + || dz , and Ly, L 2 and L 3 are Lame coefficients 

Ly = ai(l + — ), L 2 = o;2(l + -5-), L 3 = 1 (2.40) 

xiy ri 2 


2.4.1 Displacement Field 

By using the equivalent material properties, the sandwich panel may be modelled 
as a common laminated panel with the honeycomb replaced by an equivalent solid 
lamina. The geometry of the panel would be unchanged. The displacement field for 
the sandwich panel satisfying the higher order shear deformation theory is [19,20] 

u(x, y, z, t) = u(l + ^) + Z7p x + z 2 i x + z 3 Cx 

v(x,y,z,t) = vil + ^+zipy + z^ty + z*^ (2.41) 

w(x,y,z,t) — w 
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where t is time, (u,v,w) are the displacements along the (x,y,z) coordinates, (u,v,w) 
denotes the compounding displacements of the point (x,y,0) on the mid-plane and ip x 
and ip y are the rotations of normal to the mid-plane to x and y-axes respectively. 

The functions £ x ,Cx,€y and will be determined using the condition that trans- 
verse shear stresses t xz and r yz vanish on the top and bottom surfaces of the shell. 


Txz (x,y,±^,t) = 0 


T yz (z,y,±|,i) = 0 


(2.42) 


These conditions are equivalent to the requirement that the corresponding strains 
be zero on these surfaces. The transverse shear strains of a shell with two principal 
radii of curvature are given by 


€5 — 


e 4 


du 1 dw_ 
dz ' dx 

du j 1 dw 

dz ' a 1 dx 


u 

Hi 


^ + ipx + 2z£x + 3z 2 Cx + 57 

dv 1 dw_ 
dz ' dy 

dv j 1_ dw_ 

dz ' ai dy 


v 

h 2 






u 

Hi 


(2.43) 


Setting e 5 (:r, y, ±| and 64 (a;, y, ±|, to zero, we obtain 


& = 4 = 0 


Cy = ~w(^y + iS) 

Substituting eq (2.44) into eq (2.41), we obtain 


u(x,y,z) = (l + -^)u + zip x -z*^ + 
v(x,y,z) = (l + ^)v + zxp y -z 3 ^ |Vs, + ifj] 


(2.44) 


(2.45) 


This displacement field is used to compute the strain and stress, and then in 
formulating the system equations 


2.4.2 Stress-Strain Relation 

Substituting eq (2.45) into the strain-displacement relation 
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(2.46) 


e n = ei = e° + z(k ° + z 2 k\) 

€22 = €2 = ^2 -t- z(k® -t- z 2 k 2 ) 

e 33 = £3 =0 

723 = £4 = e° + z 2 k\ 

731 = e 5 = e° + z 2 k\ 

712 = e 6 = eg + + 2 2 A;|) 

where 



(2-47) 


Here x and y denote the Cartesian coordinates (dx = a x dz, dy = a 2 dz). 

The transformed stress-strain relation for the k th lamina with respect to the lam- 
inate coordinates are 



Rearrange eq (2.48) 


<2n Q 12 0 0 0 

Q21 Q22 000 

Q31 Q32 Q33 0 0 

0 0 0 Q44 q" 45 

000 q 54 q 55 
Qgi Qg2 0 0 0 




Q11 Q 12 Q i 6 
Q21 Q 22 Q26 
Q 61 Q 62 Q&G 


Cn 

^22 

712 



( ^23 ^ Q44 Qi5 

{ T 31 J {k) Qbi Q 55 


723 

731 



(2.48) 


(2.49) 


(2.50) 
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Here 


^33 — Qn^ll + Q 32^22 + Q 33633 + <536712 


au ~ ° xx a 12 -Oyy CT 33 — a zz 7"12 = T xy T 23 = T yz T 31 = r 2I 

6u = 6xx 622 = e vv e 33 = e 2Z 712 = j xy 723 = l yz 731 = 7** 

2.5 System Equation For Sandwich Panel 

The virtual variation in the total strain energy must be equal to zero. Hence 

\ ft \ rh/2 (rarb 


1 f c r ft /2 1 r a f° \ 

2 Jo J-h /2 I Jo Jo ^ Tl ^ ei + a2 ^ e 2 + a 6^6 + o' 46 €4 + o 5 6 e 5 ] dxdy — qSwdxdy > dz dt 


5 L KM l /> [(Of + (*) 2 + (*) 2 ] dxdyldz 


(2.51) 


= l [I l [ I**! + MiSKt + P 1 6K 2 + N 2 Se° + M 2 6K° + P 2 6K 2 2 + N 6 6e° 

+ MqSK® + PqSKq + Q 2^^4 + K<25K\ + Qi5e\ + K\8K\ — qSw] 
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Where 5 indicates the virtual variation in the parameter. 

Where q is the distributed transverse load, N i ,M i ,P l ,Q 2 ,Q u K 2 ,K l are the re- 
sultants, 

n rh/2 ... , 

(N h M„P,) = Y, f, oi*’(l ,z,z 2 )iz (i = 1,2,6) 

t— i J —h/2 ' 


(Qi,K 1 ) = ± [ h '\M(l,z 2 )dz 
k=i J ~ h n v 1 

(<52, ^2) = Y 1 J _ h/2 <J l fc) (1, ^ 2 ) ete 

The inertias / l andl/(i=l,2,3,4,5) are defined by the equations, 

n rh/2 

{h,h, h, h, h, h ) = X) / p (fc) (l, -z 2 , z 3 , z 4 , z 6 ) 

1 . 1 J—h/2 v J 


(2.52) 


k=1 J-h/2 
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(2.53) 


U = h + 

h = h + ~ £sh ~ uhcjs 

T^h + ih-M-ufah 

h = 3^/4 + 3 rn^h 

i* = 3p- J 4 + 3^17/5 

h = h ~ 3I2/5 + 5^/7 
I^h-ith + ^h 
h = 3^5 - ^Uz 

f ' 4 r 16 r 

J5 -3h?h~9hi-l7 

The governing equations of motion can be derived from eq (2.51) by integrating the 
displacement gradients by parts and setting the coefficients of 5u, 6v, Sw, 8ip x and 8ip y 
to zero separately with the moment terms in the first two equations being omitted 
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(2.54) 

r 8w 
dx 
,dw 

dy 


The Resultants are expressed in terms of the strain components using eq(2.46) 


and eq(2.48)in eq(2.52) we get 


N t = + B tJ K° + E^K 2 

M t = B l3 e] + D tj K J + FijK 2 (*, j = 1, 2, 6) 

Pi = E i3 e° 3 + FijK J + HijKj 

Q 2 = A 4j e° 3 +D 4j K} (2.55) 

K 2 = D 4j e 0 j+F 4j K j (j = 4, 5) 

K\ — D 5j e°j + FsjKj 
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Where A t j, B t j, 1 

(A t j, 

for i,j=l,2,4,5,6. 

dNi _ 

dx 


dN 2 

dy 

dNe 

dx 

dNe 

dy 

dMi 

dx 

dM 2 

dy 

dM 6 

dx 

dMe 

dy 


Aj j Eij 5 Bij , H tJ are the laminate stiffnesses 
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For a curved panel having a rectangular plane from, the Simply supported boundary 
conditions may be expressed as 
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at x=0,a 


v(0,y) = v(a ,y) = 0 
w(0,y) = w(a,y) = 0 
N i(0,y) = Nx(a, y) = 0 
(0,y) = Mi (a, y) = 0 
A(0 5 y) = Pi(a,y) = 0 

Vv(°»y) = Vv(°>y) = o 

at y=0,b 


(2.58) 


u(x, 0) = u(x,b) = 0 
w(x, 0) = w(x, b) = 0 
N 2 (x, 0) = N 2 (x, b) = 0 
M 2 (x, 0) = M 2 (x, b) = 0 
Pi(x,0) = P 2 ( x ,b) = 0 

ipx(x,Q) = tp x (x,b) = 0 

Where a and b denote the length along the x-and y-direction, respectively. For a full 
cylinder the boundary conditions on X are not there. 

We assume the following Naviour type solution form that satisfies the boundary 
condition in eq (2.58) 


/ ,\ x — \ v — r TT TTITTX . TVKTJ 

u {x,y,t) — 2 ^ > , U mn cos sm-~- sinuo mn t 

m= 1 n — 1 ^ 0 

oo oo 

/ . TflTT X Tb'KXj 

v \X,y,t) — 2_s / , VmnSin cos—M-sinui mn t 

m = 1 n=l a b 

OO OO 

( j.\ X~^ V "' 1 t xt * TTITTX . TITTV 

w[x, y,t) = ^ 2_, w mn sm S zn~-sinoj mn t ( 2 . 59 ) 

m=in=i a b ' 

OO 00 

m~l n=l ^ 0 

OO 00 

= E E &.«»— 

m=l n=l G b 

Now we Fourier Transform the displacement function split in the three parts 

/ +°° rO /*1 r-l-oo 

/(*)=/ m+fm+f m 

-oo •/—oo JO 7 1 

Hsrc integral oo to 0 is zero because the displacement is undefined in the negative 
domain of the set of integers. The integral 0 to 1 is zero because m and n are not 
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defined for the range 0 to 1. Where N is maximum time limit, Substituting eq(2.59) 
into eq(2.54) we obtain 


' M n 

M i2 

CO 

m 14 

m 15 


r .. \ 

U mn 


^11 

&12 

^13 

&14 

kis 


^mn 


M 2i 

M 22 

M 23 

M 24 

M 25 


Knn 
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&22 

&23 

&24 

k 25 


hmn 


CO 

M 32 

M 33 

M 34 

M 35 

< 

Wmn 

► + 

^31 

&32 

&33 

&34 

& 35 

< 

W mn 

> == < 

M n 

M 42 

M 43 

m 44 

m 45 


^ mn 


kn 

^42 

&43 

^44 

& 45 
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_m 51 

M 52 

M 53 

M, 54 

M 55 _ 


k r mn 
> > 


&51 

&52 

&53 

k $4 

& 55 


ih'k 

. > mn 

\ j 



(2.60) 

for any m , n 

Where Q mn are the coefficients in the double Fourier expansion of the transverse 
load, 

OO OO 

„ m 7TX . TiTTV . 

-smuj mn t (2.61) 


= E YLQrnnSin 

m=l n= 1 


- sin - 


a b 

and the coefficients M tJ and k io (*, j = 1, 2 , . . . , 5) are given in the Appendix 

Equation(2.59) can be solved for U mn , V mn , W mn , ^ n , for each m and n, and 
then the solution is given by eq(2.58) are evaluated using a finite number of terms 
in the series. For free vibration analysis, eq(2.59) can be expressed as an eigenvalue 
equation 
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k 32 

&33 
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^ mn 
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m 33 

m 34 

M 35 
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w 
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k 4i 
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^43 

^44 

k 45 


m 41 
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m 43 

m 44 
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k 52 
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_ M 51 

M 52 

M 5 3 

m 54 

M55 J 
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(2.62) 


0 

0 

Qmn 

0 

0 


Where oo is the frequency of natural vibration 


2.6 Composite Material 

From eq(2.49-2.50) re write Stiffness Matrix 




Qll Q 12 Ql6 


' \ 

^11 

^22 

► = 

Q 21 Q 22 Q26 

< 

e 2 2 * 

7*12 

V > 

(*) 

Q 61 $62 Q66 

(*) 

T 12 

V. / 


35 


(2.63) 


(fc) 



(2.64) 


U] 


Q44 Q45 


f 723 

< 

1 



< 


l r 31 J 


Q54 Q55 

(*:) 

{ 731 J 


Here 

Ql6 = $61 
Q26 = Q&2 
Q45 = Q54 

Qn,Q 22 , Q i6) Q26, Qm, Q 44 , and (Q 5S we can get from Ref [21]. Substituting eq (2.62 
to 2.64) in eq (2.51) after simplifying we can get response often satisfying boundary 
condition. 
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Chapter 3 


Solution Approach 

3.1 Eigen Value Problem 

The System is assumed to vibrate in the principal mode with the displacements as ex- 
pressed in eq (2.59). where the shape functions u(x, y, t),v(x, y, t), w(x, y, t ), ip x (: r, y, t ) 
and ip y (x,y,t ) satisfy the boundary conditions for the shell. 

The Shape functions are initially assumed to be zero then constituted to wave 
forms. These can be expressed as in eq (2.59). where U, V, W, ip 1 and ip 2 are the dis- 
placement amplitudes sinusoidal functions with (^f 2 ) and ( 2 p) as their arguments, 
respectively. Here m and n denote longitudinal and circumference wave numbers. 

Eq (2.62) can be written in matrix notation 


([if] - w 2 [M]) {A} = 0 

(3.1) 

Where K and M are (5 x 5) square matrices and 


A = V mn Wmn ^mn t^Vnn 

T 

(3.2) 

Here 

[M] t = [M] cl + [M] h + [M\ c2 


[K] t = [K] cl + [K] H + {K} c 2 

(3.3) 

And 

[P] = [M]^[K] t 

(3.4) 
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The subscription “c” denotes laminated panel and “H” denotes the honeycomb panel. 
Subscription “t” denotes total mass and stiffness matrices of the panel. Equation (3.1) 
can be expressed as 

([P] - A[i]) {A} = 0 (3.5) 

Where {A} is the amplitude vector and A = co 2 . The eigen value problem becomes 

|[P]-A[/]|=0 (3.6) 

Expanding eq (3.6) we get a fifth order characteristic equation 

A 5 - B x A 4 - B 2 A 3 - B 3 A 2 -B 4 X-B 5 = 0 (3.7) 

Where B X ,B 2 ,B 3 , B 4 and B 5 are the coefficient and are placed in Appendix [B]. 
Prom eq (3.7) for the deterministic problem eigen value solution is obtained in the 
usual manner with standard technique. However, the eigen value A as well as the 
eigen vector {A} are random as the characteristic equation has random coefficients. 
The eigen value and eigen vector statistics in this case can not be found directly and 
more analysis is needed. A procedure is being outlined here. 

3.2 Perturbation Approach 

We write the random components as the sum of the mean value and the zero mean 
random part. 

A„ = A™ + A;; {&„} = { A™} + {A;} ; [K] = [*]"> + [K] r n = 1, .... 5 (3.8) 

It can be seen that 

A™ = «) 2 and A; = 2wJX + «) 2 n = l,...,5 (3.9) 

The subscription “r” denotes zero mean random part and “m” denotes the mean part 
of the parameters. 

It is assumed that the random component are very small compared to the mean 
value for each parameter. Substituting eq (3.8) and eq (3.9) in eq(3.1) and collecting 
the terms of zero and first order small quantities 

[ir] m {A”} = A"[M]{A”} 
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(3.10) 



and «J£P - A“[M]) {A;} = - ([*]' - A;[M)) {A™} (3.11) 

Where [K] and [M] denote the total system quantities. Eq (3.10) is in the mean and 
the mean values of the eigen value and eigen vectors axe obtained in the usual manner 
with the help of any standard procedure. 

The eigen vector are orthogonal with each other and would satisfy the condition 

{A”} T M{Ar}=^ and {AJ>} r [K] {A|»} = A™<5„, n = 1 5(3.12) 

All eigen values are assumed to be distinct. In the given space the eigen vectors 
form a complete orthogonal set and any vector in the space can be expressed in terms 
of these. Hence, the random component 

{A;}=Eftu{A,} m Cnn = 0 (3.13) 

1=1 

Here c n i are unknown random coefficients and have to be evaluated. 

Substituting eq (3.13) in eq (3.11) 

£ (M m - a“[m]) c* {ah = - ([*r - a;[m]) {a™} ( 3 . 14 ) 

l - 1 

Premultiply by {A", eq (3.14) becomes 

£ {A”} T ([K] m - A ”[M)) {AH = {A"f (IK]' - AJ[M]) {A"} 

i = 1 

or £ [{A“} T [KT {A?*} C,a - {AJ} T [M] {AD] 

1 = 1 

= + {A"f [. K ]' {A"} - A; {A-} t [M] {A?} 

or £ [A - A„ Cnidni] = + { A“} T [K] r {An } - A;i 

1 = 1 

or A; = {A™f IK]' {A™} (3.15) 

Premultiplying eq (3.14) by {A”} r with k^n 

£ {A?f ([X]™ - A ”[«]) Cni {AD = {A?f ([KY - A„[M]) {AJ} 

1=1 

£ [{ Ai“f [KT {AD OW - A^c„i { A?} t [M] {Ar}] 

1 = 1 

= + {A£>} r m {An - a; {AH r [M] {A-} 
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or 



“ r E WMm<W - = + {A? } r [K\ r {A”} - Wb, 

1=1 


°r M (A? - A”) = + {A?} r IKY {A“} 


, , _ {AJ*} [KY {A™} 

L J - \™ 


(3.16) 

using eq (3.16) in eq (3.13) 



{AU = E {Ar ^ M (i A " } (AD 

1=1 A l “ A n 

l 77 1 

(3.17) 


Variance of A n and {A n } can be obtained from eq (3.15) and eq (3.17). Consider 
the RHS of eq (3.15) and let 


{ A n } [QiiI Qn2 <5n3 Qn4 Qn§\ 

Then 

{A”f [KY {a™} = EE qZkic, = E E «3«3*& 

1=1 j=i 1=1 i=i 

Hence from eq (3.15) 

{a ;} 2 = E E iXfij E E 

(=1 J= 1 S =1 t=l 

«• R } 2 = E E E E tXjCCMK, 

1=1 3=1 3=1 t=l 

Var {A;} = EEEE 9 SJlK" (*&*-) (3.18) 

1=1 J =1 3=1 t=l 

The zero mean random parts of elements of stiffness matrix are expanded in Taylor 
Series. Using a single step Taylor series in-terms of the basic material properties likes 
-Eii, #22 etc, 


k r 

K ij 


v dk ? 1 v fll?r 

Y' OK v at i y —ii A r 

z+ ^dA, m ‘ 1 


(3.19) 


Where Ai are the basic Material properties and rj is the number of basic material 
properties. The second and higher order terms in eq (3.19) has been dropped because 
there are very small compared to the first order terms. Further, the second set of 
terms are very small compared to the first set and may be ignored in evaluation. 
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3.3 FEM Approach 

The random problem can be solved by using finite element method (FEM) with Monte 
Carlo simulation. In this approach the structure is discretised in to finite elements 
with nodes. The system properties and external loading at each node is obtained by 
computer simulation of numbers that first this given statistical distributions. The 
problem in solved repeatedly to obtain a sample of the response statistics. This 
sample is analysed to obtain the response. 

These steps are being described below for the present problem 

3.3.1 Generating Random Number 

The material properties, used as basic variables, have been generated assuming them 
to be Gaussian. The random numbers are generated first in the range -1 to +1 in the 
normal distribution curve. After that for analysis we are assuming the basic material 
properties to have a input mean value and input of standard deviation value. The 
generated numbers are modified to have the desired mean and SD. New basic material 
properties are generated for each cycle of solution. These basic material values are 
the input value of the FEM package for repeated solution. 

3.3.2 Modelling 

Geometric modelling has been performed using MSC-Patran package. The geometric 
model used is shown in Fig (3.1). This model has three set of surfaces. The top 
and bottom surfaces are laminated surfaces while the middle surface is honeycomb 
structure. 

The most distinctive feature of the finite element method that separates it from 
others is the division of a given domain into a set of simple sub-domains called finite 
element. Any geometric shape that allows computation of the solution or approxi- 
mation can be used as a finite element. Finite element discretization is represented 
as a collection of a finite number of n sub-domains, namely line, surface and volume 
segments. Each sub-domain is called an element. The collection of elements is called 
finite element mesh. The elements are connected to each other at points called nodes. 
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FEM mesh used for this problem contains approximately 4600 sub-domain each 
sub-domain contains four noded quadratic rectangular shell elements. After generat- 
ing elements in different parts of the structure it is required to merge them together. 
This becomes necessary as this model contains different surface and different material 
properties. When we are generating elements it will created two node identity in the 
same location. For avoiding this we need to do an equivalence operation in the finite 
element model. After merging the finite element model it becomes a single geometry. 
Each nodes has been assumed to possess 5 degrees of freedom-three translational and 
two rotational. The FEM model is shown in Fig (3.2) 
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Figure 3.2: FEM Model of Honeycomb structure 
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Chapter 4 


Results And Discussion 


The solution approaches discussed in the previous chapter have been used to obtain 
the natural frequencies of honeycomb sandwich laminated cylinder, simply supported 
at the two ends. The sandwich is assumed to have facing laminates on the two sides 
with a hexagonal honeycomb core layer. The facing laminates have each lamina of the 
same thickness. The core is of aluminium and the facing is taken to be graphite epoxy. 
The mean material properties adopted for generating the results are Material proper- 
ties for equivalent solid aluminium core: E a = 72395 N/mm 2 , G s = 27578.8 N/mm 2 , 
= 0.3125, p = 2.79567 * 10-6 Kg/mm 3 . 

Material properties for Graphite/epoxy^n = 134000 N/mm 2 , E 22 = 10200 N/mm 2 , 
G w = 5520 N/mm 2 , G 23 = 3430 N/mm 2 , G n = 5520 N/mm 2 , v 12 = 0.3, 
p = 1.577 * 10~ 6 Kg/mm 3 . 

The natural frequencies have been nondimensionalized as 


<N I 

cdl 

/ p 

M 

1 E22 


Where u is the natural frequency 

R is the radius of the cylinder. 

h is the total thickness of the laminate 

p is the density of Graphite/epoxy 

E 22 is the material properties of Graphite/epoxy 

u is the nondimenalized natural frequency 
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4.1 Validation 


Validation has been sought for the approach by comparison of results with published 
data. For this, only the mean value has been compared as information on variance is 
not available in literature. Table (4.1) contains non dimensionalized natural frequency 
obtained by the present work and its comparison with Ref [20]. It can be seen that 
the results exhibit a close agreement. 


4.2 Natural frequency statistics 

The perturbation approach discussed in the chapter 3 has been used to find the mean 
and the variance of the first five natural frequencies. Three stacking sequences for the 
facing layer: (l)-[0790°/90707iy/0790 o /90 o /90°], (2)-[0790707907ff/0790790790°] 
and (3)-[90 o /0°/iii r /0 o /90 o ] have been studied for different length to radius (L/R) ra- 
tios and different thicknesses. Behaviour of the nondimensionalized mean and vari- 
ance are presented below. 

4.2.1 Mean frequency 

Table (4.2 to 4.7) present the first five natural frequencies for different length ratios. 
There are different lay up sequences with each layer thickness 0.1 mm. Mode shape 
number m and n are selected from 1,2 and 3 in all possible combination. When we 
compare first mode shape m=n=l for different L/R. We see that non-dimensional 
natural frequencies decrease as L/R increases. This may be because stiffness decreases 
and mass increases with increases in L/R. Comparing mode shapes m=l,n=2 and 
m=2,n=l it can be seen that for m=l,n=2 non dimensional natural frequency is 
higher. It indicates that change in mode shape in y-direction is dominant for natural 
frequency. This behaviour is common to men other layup combinations. 

The natural frequency is seen to increase with thickness, as expected, for all the 
three layup sequences. The third layup sequences has the highest natural frequency 
while the second has the smallest. The anti-symmetric case has [B l0 ] ^ 0 and it may 
effect this region for increase or decrease in non dimensional frequencies. Table 4.4 
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also presents a symmetric lay up sequence case. But the lay up has lower number 
of layers compared with the other tables. The general behaviour discussed above are 
the same for the Tables (4.5) to (4.7) where the honey comb thickness changes but 
the L/R is constant. 

4.3 Natural Frequency Variance 

The results are presented graphically for changes in the square of SD/mean of the 
natural frequencies with changes in SD/mean of the input variables. Study has been 
conducted for the three layup sequences. However the results are presented only for 
the first layup in details because of the high space requirement. The study has been 
conducted first with each input variable changing alone, one by one, and then with 
all of them changing simultaneously. The variation in input SD/mean is taken from 
0 to 20%. 

Fig (4.1) to (4.27) for the first layup and figure 4.28 for the second layup show that 
the variance of natural frequencies change linearly with the basic material property 
SD. The main reason for this may be due to the inherent linearsation in the first order 
perturbation approach as well as in the one term Taylor series expansion employed 
in the analysis. 

4.3.1 Only Changing 

Fig (4.1) to (4.3) shows the effect of randomness in En on the variance of the natural 
frequency. E r n has a most pre dominant effect on the SD of the square natural 
frequencies as compared to the other input variable. More over as seen from these 
curves first natural frequencies are dominant compared to the 2 nd and higher natural 
frequencies. This indicate that the effect of randomness is very high in the first mode. 
When we consider increasing mode shape natural frequency variances are reducing 
compared to the first mode shape. Comparison of m=l,n=3 and m=3,n=l shows 
that y-direction mode shape has dominant effect on natural frequency variance. It is 
indicated that effect of randomness is very high in that direction. With L/R increase 
non dimensional SD squared natural frequencies increase. The general behaviour for 
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the lay up sequences 0°/90°/0 c, /90 o /i?'/0 o /90 o /90 o /0 o is similar except that there is 
differences in frequencies. 

4.3.2 Only E 22 Changing 

Fig (4.4) to (4.6) show that the variance of the natural frequencies increases linearly 
with the input of random value E 22 . Comparison indicates that for mode shape 
m=l,n=l the first five non dimensionalized natural frequencies variances are domi- 
nant. Other natural frequencies have very low variances compared to the first nat- 
ural frequencies. Another combination of mode shape m=l,n=l,2,3 the first mode 
shape m=l,n=l also indicates a higher variation. It’s valid through out the com- 
bination of mode shape and L/R ratio. The behaviour for the lay up sequences 
0 0 /90°/0°/90 0 / jy r /0°/90 o /90°/ 0° is also similar to the first lay up except for differ- 
ences in frequency values. 

4.3.3 Only G 12 Changing 

Fig (4.7) to (4.9) show that variance of the natural frequencies increases linearly 
with the SD of input random value Gi 2 . First natural frequencies variance is again 
dominant compared to second and higher natural frequencies. This observation is 
valid for different combination of the mode shape and L/R ratio. The dynamic 
behaviour for the lay up sequences 0°/90 o /0 o /90 o /fl’/0 o /90 o /90 o /0 o is also similar. 

4.3.4 Only G 23 Changing 

Fig (4.10) to (4.12) show that the variance of the natural frequencies increases lin- 
early with the input random value G 2 3- Details follow the previous results. For the 
4 th and 5 th natural frequencies values are very close for all mode shapes. This in- 
dicates that shear effect due to G 23 is low. This is valid for different combination 
of mode shapes and L/R ratios. The dynamic behaviour for the lay up sequences 
0 0 /90 0 /0 0 /90 0 /JT/0 0 /90 0 /90°/0 0 is again similar. 
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4.3.5 Only G 31 Changing 

Fig (4.13) to (4.15) present the natural frequency variance with G 31 only changing. 
The results follow the pattern of G 2 3 variation. Again the dynamic behaviour for 
the lay up sequences 0 0 /90°/0 0 /90°/iif/0 0 /90 0 /90 0 /0 0 is similar with only differences 
frequencies having different values. 

4.3.6 Only E s Changing 

Fig (4.16) to (4.18) variance are for E s only random. The natural frequencies of 
structure mode shape m=l,n=3 are the dominant natural frequencies. The effect 
is very high in the y-direction. The second and third non dimensionalized natural 
frequencies are above of the first natural frequency. These are for LJR = 1.5. The 
nature of the graph change with change in L/R ratio. 

4.3.7 Only G s Changing 

Fig (4.19) to (4.21) present the results for only G 3 having randomness. The fun- 
damental natural frequencies of honeycomb structure mode shape m=l,n=l is the 
dominant natural frequency. The effect is very high in the y-direction. The above 
behaviour changes with decrease in the L/R ratio. The frequency variance for the 
lay up sequences 0 o /90 o /0 o /90 o /ff/0 o /90 o /90 o /0° exhibits a similar nature with only 
differences in magnitude of the frequencies. 

4.3.8 Only z /12 Changing 

Fig (4.22) to (4.24) shows the results for ^ 12 only varying. The fundamental natural 
frequency variance for mode shape m=l,n=l is the dominated natural frequency. The 
effect is very high in the y-direction. For the second and third natural frequencies the 
order of sequence of magnitudes of the mode shape is reversed with respect to the 
first natural frequency. With decrease in the L/R ratio the dominant mode changes, 
m=2,n=l being dominant for L/R= 1.0 and m=3,n=l for L/R — 0.5. 
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Table 4.1: Comparison of mean u 


Lay Up sequence 

Present work 

J.N. Reddy [20] 

a/h=10 

a/h=100 

a/h=10 

a/h=100 

0°/90° 

8.978141 

9.712483 

8.975 

9.712 

0°/9070° 

11.71697 

15.039365 

11.79 

15.19 

0°/90°/90°/0° 

11.776354 

15.188267 

11.78 

15.19 


4.3.9 All Material Properties Simultaneously Random 

The variance of the natural frequencies are plotted against the SD of material prop- 
erties in Fig (4.25) to (4.27). The compared fundamental natural frequencies of 
honeycomb structure mode shape m=l,n=l is the dominant natural frequency. The 
y direction mode has very high effect on the natural frequency variance. If compared 
to the first second and third non dimensionalized natural frequencies have reverse 
order of magnitudes. With change in L/R ratio the above behaviour changes for the 
honeycomb structure is shown in Fig (4.26) and (4.27). The dynamic behaviour for 
the lay up sequences 0 0 /90°/0 0 /90°/iif/0 0 /90 0 /90 0 /0 0 is similar. 

4.4 FEM Result 

Fig (4.28 ) shows the plot using the FEM with Monte-Carlo simulation. A sample 
sets of 20 was analysed for the result, the large size of the problem requires 4 hrs 
for one point. The sample size, however is not large enough to be representative as 
indicated by the nature of the case and it was not possible to generate a large sample. 


k~T f'-TT-T #=7 TIT 


.. . rv.. r r — 

s»io A- 
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Table 4.2: Mean natural frequency for different L/R , lay up sequence(l) 


Mode Shape 

0°/ 90°/90°/ 0°/ H/0° /90° /90° /0° 


L/E=1.5 

L/R= 1.0 

L/R=0.5 


w„i = 

uJ t j i = 1 , ... 3 5 

u t ,i = 1, — ,5 

m=l,n=l 

227.171005 

327.643420 

550.343644 

1251.699907 

1620.863005 

2173.776609 

2068.663326 

2485.669475 

4491.111449 

223180.769586 

223180.780103 

223180.836911 

231671.576074 

231680.442009 

231728.307533 

m=l,n =2 

415.111843 

550.344016 

767.618573 

1828.939186 

2173.777820 

3421.960019 

3101.208623 

4491.111425 

8836.083742 

223180.794828 

223180.836911 

223181.064470 

231692.853281 

231728.307517 

231919.595005 

m=l,n=3 

550.344385 

680.455598 

906.860374 

2173.779033 

2760.210425 

4856.601120 

4491.111399 

6653.311944 

13218.327575 

223180.836911 

223180.931665 

223181.444818 

231728.307501 

231808.043989 

232237.777910 

m= 2 ,n=l 

358.409705 

528.029733 

954.138765 

1334.313420 

1791.958327 

2374.885245 

2486.778188 

2760.141615 

4565.429197 

223182.286270 

223182.295178 

223182.343359 

231671.578847 

231680.446282 

231728.319845 
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Table 4.2: (Continued) 


Mode Shape 

0°/90 o /90 o /0 o /H/0 o /90 o /9Q° /0° 


L/R=1.5 

L/R=1.0 

L/R=0.5 


w„i = 

u) x 3 i 1 j . . . j 5 

u> t ,i = 1, — ,5 


686.074178 

954.139613 

1369.390039 


2062.175534 

2374.886212 

3476.465713 

m=2,n=2 

3256.809887 

4565.429114 

8866.241038 


223182.307656 

223182.343359 

223182.537385 


231692.859647 

231728.319829 

231919.638507 


954.140461 

1221.088877 

1527.294495 


2374.887182 

2867.108582 

4874.474918 

m=2,n=3 

4565.429029 

6695.709334 

13237.700387 


223182.343359 

223182.423954 

223182.865093 


231728.319813 

231808.069482 

232237.870168 


423.728014 

631.155806 

1208.809303 


1432.008358 

1948.058317 

2721.534733 

ra=3,n=l 

3085.675040 

3251.726421 

4710.048147 


223184.814035 

223184.820257 

223184.854056 


231671.583469 

231680.453395 

231728.340373 


832.887208 

1208.810570 

1937.991624 


2337.101310 

2721.535494 

3596.809944 

m=3,n=2 

3582.331018 

4710.047961 

8918.417292 


223184.828993 

223184.854056 

223184.992168 


231692.870262 

231728.340357 

231919.710892 


1208.811834 

1659.094003 

2185.491181 


2721.536254 

3098.280886 

4909.992944 

m=3,n=3 

4710.047776 

6771.353630 

13270.516835 


223184.854056 

223184.911043 

223185.232117 


231728.340342 

231808.111989 

232238.023992 
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Table 4.3: Mean natural frequency for different L/R, lay up sequence(2) 


Mode Shape 

0 0 /90 0 /0 0 /90 0 /H/0 0 /9Q 0 /90°/Q 0 


L/R=1.5 

L/R=1.0 

L/R=0.5 


uJj,i= 1, — ,5 

Ui,i = 1,...,5 

Wf,i = l,...,5 


222.855516 

321.977084 

544.304261 


1250.818722 

1616.575003 

2153.042833 

m=l,n=l 

2066.958436 

2484.137519 

4490.636360 


225332.012723 

225332.022608 

225332.076060 


229579.785761 

229588.760073 

229637.210374 


408.790928 

544.304635 

764.832394 


1819.690241 

2153.044004 

3364.004657 

3 

II 

l—i 

11 

to 

3100.275615 

4490.636335 

8835.886370 


225332.036453 

225332.076060 

225332.291109 


229601.323027 

229637.210359 

229830.829400 


544.305006 

676.197275 

905.595159 


2153.045178 

2720.699332 

4763.892141 

CO 

11 

t-H 

ii 

a 

4490.636308 

6653.035605 

13218.200537 


225332.076060 

225332.165430 

225332.653400 


229637.210343 

229717.919771 

230152.871929 


350.967391 

517.506946 

939.571039 


1334.012192 

1790.236093 

2360.488803 

m=2,n=l 

2485.742301 

2758.684607 

4564.896553 


225333.526286 

225333.532981 

225333.569442 


229579.789819 

229588.767210 

229637.233906 
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Table 4.3: (Continued) 


Mode Shape 

0°/9070790 o /if/0° /90° /90° /0° 


L/R=1.5 

L/R=1.0 

L/R=0.5 


j % — 1 j ■ - . j 5 

uJi : i = 1 , . . . , 5 

= 1, ... ,5 


673.235105 

939.571886 

1361.763647 


2057.264838 

2360.489735 

3421.869080 

m=2,n=2 

3255.720707 

4564.896468 

8866.037444 


225333.542383 

225333.569442 

225333.720244 


229601.334453 

229637.233891 

229830.914802 


939.572732 

1209.495359 

1523.886919 


2360.490670 

2833.082808 

4782.964179 

m=2,n=3 

4564.896382 

6695.417502 

13237.571576 


225333.569442 

225333.631341 

225333.987106 


229637.233876 

229717.969781 

230153.048823 


414.693250 

617.911182 

1186.211460 


1431.911412 

1947.530370 

2714.059578 

m=3,n=l 

3085.117285 

3250.700391 

4709.419210 


225336.048845 

225336.050216 

225336.058337 


229579.796586 

229588.779104 

229637.273160 


815.842923 

1186.212711 

1921.600647 


2335.246843 

2714.060317 

3550.073743 

m=3,n=2 

3581.204860 

4709.419023 

8918.202909 


225336.052210 

225336.058337 

225336.101976 


229601.353510 

229637.273144 

229831.057103 


1186.213960 

1636.048599 

2178.819924 


2714.061057 

3074.656795 

4821.105997 

m=3,n=3 

4709.418836 

6771.034502 

13270.385004 


225336.058337 

225336.074415 

225336.209699 


229637.273129 

229718.053203 

230153.343856 
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Table 4.4: Mean natural frequency for different L/R, lay up sequence(3) 


Mode Shape 

90 o /07 H/0°/90° 


L/R=1.5 

L/R=1.0 

L/R=0.5 


cjf j i = 1 , . . . , 5 

u t ,i = 1,. . . ,5 

w„i = 1,...,5 

m=l,n=l 

277.546514 

400.447120 

717.440270 

1231.512680 

1469.653621 

1763.137796 

1983.809895 

2493.473141 

4600.358711 

215540.014275 

215540.061561 

215540.146558 

215548.049243 

215558.295447 

215613.779456 

m=l,n=2 

517.546930 

717.440813 

1036.250041 

1580.321471 

1763.138499 

2586.442502 

3159.837519 

4600.358686 

9059.410681 

215540.090374 

215540.146558 

215540.408666 

215572.675968 

215613.779447 

215835.529846 

m=l,n=3 

717.441356 

918.549642 

1182.996255 

1763.139197 

2124.539398 

3656.147965 

4600.358663 

6820.772003 

13552.687539 

215540.146558 

215540.257431 

215540.837931 

215613.779438 

215706.224359 

216204.177258 

m=2,n=l 

388.709035 

536.205444 

1012.643294 

1496.840694 

1983.914046 

2436.538841 

2718.580112 

2927.011548 

4698.708637 

215539.896063 

215540.015533 

215540.156070 

215549.912872 

215560.086736 

215615.514365 
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Table 4.4: (Continued) 


Mode Shape 

m o /0°/H/0°/90 o 


L/R=1.5 


L/R=0.5 


w„* = 1,.. .,5 


uJj, z = 1, . . . j 5 


694.780152 

1012.644323 

1787.200793 


2254.472189 

2436.539218 

2848.534304 

m=2,n=2 

3381.967610 

4698.708542 

9098.474611 


215540.075564 

215540.156070 

215540.433604 


215574.435821 

215615.514356 

215837.245979 


1012.645350 

1450.989842 

2108.551425 


2436.539593 

2598.232382 

3727.795131 

m=2,n=3 

4698.708448 

6875.885070 

13577.739173 


215540.156070 

215540.278307 

215540.865816 


215615.514346 

215707.946505 

216205.884905 


469.595400 

607.656444 

1109.196051 


1678.407943 

2245.106409 

3260.350769 

m=3,n=l 

3719.319095 

3803.910107 

4926.335054 


215539.773959 

215539.955241 

215540.171153 


215552.943949 

215563.055726 

215618.406593 


768.248933 

1109.197199 

2128.802737 


2749.994214 

3260.351197 

3549.191687 

3 

II 

CO 

'b 

II 

bO 

3992.383691 

4926.334819 

9167.988489 


215540.054336 

215540.171154 

215540.474576 


215577.365398 

215618.406584 

215840.106547 


1109.198346 

1627.871529 

2905.284143 


3260.351621 

3442.606534 

3978.288363 

m=3,n=3 

4926.334587 

6980.108038 

13620.647249 


215540.171154 

215540.312241 

215540.911994 


215618.406574 

215710.817553 

216208.731225 
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Table 4.5: Mean natural frequency for different Web Thickness, lay up sequence (1) 


Mode Shape 

0790 o /90707if/0°/90 o /90° /0° 


Thickness=0.3 

Thickness=0.4 

Thickness=0.5 


O l} i = 1, . . . , 5 

uM = 1,.. .,5 

C0i,i = 1,. . . ,5 


227.171005 

244.069519 

263.030256 


1251.699907 

1341.239593 

1397.811743 

r—5 

II 

1 l 

II 

a 

2068.663326 

2234.064423 

2404.081059 


223180.769586 

249709.882679 

273679.367457 


231671.576074 

257328.744608 

280650.048367 


415.111843 

447.832763 

487.659151 


1828.939186 

1861.163378 

1883.327867 

xn=l,n=2 

3101.208623 

3558.285539 

3966.609262 


223180.794828 

249709.936861 

273679.464754 


231692.853281 

257354.652822 

280680.049593 


550.344385 

599.477407 

663.986306 


2173.779033 

2217.895099 

2249.234981 

m=l,n=2 

4491.111399 

5181.101772 

5792.143978 


223180.836911 

249710.027194 

273679.626972 


231728.307501 

257397.822596 

280730.038470 


358.409705 

375.622287 

390.300630 


1334.313420 

1445.739427 

1524.084937 

m=2,n=l 

2486.778188 

2633.911490 

2787.811788 


223182.286270 

249711.641153 

273681.340420 


231671.578847 

257328.749732 

280650.056759 
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Table 4.5: (Continued) 


Mode Shape 

079079070 o /i?/0 o /90° / 9070 0 


Thickness=0.3 

Thickness=0.4 

Thickness=0.5 


cj Z 5 i = 1, , 5 

cJj, i = 1, . . . , 5 

uj t ,i = 1, . . - , 5 


686.074178 

718.316880 

747.768978 


2062.175534 

2105.884091 

2133.898486 

m=2,n=2 

3256.809887 

3700.600092 

4104.941705 


223182.307656 

249711.689288 

273681.429168 


231692.859647 

257354.663656 

280680.066133 


954.140461 

997.838583 

1042.676006 


2374.887182 

2412.933294 

2441.842038 

m=2,n=3 

4565.429029 

5255.829465 

5868.882707 


223182.343359 

249711.769637 

273681.577302 


231728.319813 

257397.842856 

280730.068423 


423.728014 

442.097286 

455.942889 


1432.008358 

1554.646191 

1647.996454 

B 

II 

CO 

II 

1— > 

3085.675040 

3228.271912 

3375.990203 


223184.814035 

249714.571891 

273684.628631 


231671.583469 

257328.758273 

280650.070750 


832.887208 

868.448796 

896.048728 


72337.101310 

2423.188230 

2473.730565 

m=3,n=2 

3582.331018 

3985.552629 

4372.091336 


223184.828993 

249714.609942 

273684.703121 


231692.870262 

257354.681722 

280680.093712 


1208.811834 

1258.712205 

1299.812094 


2721.536254 

2759.720089 

2788.918189 

m=3,n=3 

4710.047776 

5395.089821 

6008.366973 


223184.854056 

249714.673641 

273684.827767 


231728.340342 

257397.876641 

280730.118375 
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Table 4.6: Mean natural frequency for different Web Thickness, lay up sequence (2) 


Mode Shape 

0790707907 ^/ 0 ° /90° /90° /0° 

Thickness=0.3 

Thickness=0.4 

Thickness=0.5 

cjj , i 1 j . . . 3 5 

u) % 1 j . , . j 5 

= 1,...,5 

B 

II 

h- 1 

'b 

II 

I— 1 

222.855516 

239.436150 

258.314104 

1250.818722 

1339.797536 

1395.895394 

2066.958436 

2232.514121 

2402.719802 

225332.012723 

251634.428104 

275436.482489 

229579.785761 

255447.144531 

278925.818106 

m=l,n=2 

408.790928 

441.217378 

481.038511 

1819.690241 

1850.152147 

1871.111758 

3100.275615 

3557.537933 

3965.968250 

225332.036453 

251634.479884 

275436.576363 

229601.323027 

255473.314661 

278956.077975 

3 

II 

h- 1 

P 

II 

CO 

544.305006 

593.369312 

658.021488 

2153.045178 

2193.992466 

2223.136989 

4490.636308 

5180.684953 

5791.767812 

225332.076060 

251634.566304 

275436.733034 

229637.210343 

255516.920544 

279006.497418 

m=2,n=l 

350.967391 

367.385401 

381.582349 

1334.012192 

1445.228674 

1523.370657 

2485.742301 

2632.872825 

2786.827774 

225333.526286 

251636.182433 

275438.450353 

229579.789819 

255447.151679 

278925.829369 
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Table 4.6: (Continued) 


Mode Shape 

0 o /90707907if/0 o /90° /90° /o° 


Thickness=0.3 

Thickness=0.4 

Thickness=0.5 


u>i,i = 1, ...,5 

% — 1 j j 5 

w„i = l,... ,5 


673.235105 

704.271626 

732.991570 


2057.264838 

2099.717910 

2126.902282 

3 

II 

to 

» 

II 

to 

3255.720707 

3699.763764 

4104.246584 


225333.542383 

251636.222244 

275438.527317 


229601.334453 

255473.333437 

278956.105742 


939.572732 

982.250223 

1026.459805 


2360.490670 

2395.986015 

2423.182690 

m=2,n=3 

4564.896382 

5255.379548 

5868.485368 


225333.569442 

251636.289081 

275438.656447 


229637.233876 

255516.958338 

279006.552050 


414.693250 

432.027639 

445.191189 


1431.911412 

1554.481064 

1647.760379 

m=3,n=l 

3085.117285 

3227.681503 

3375.399949 


225336.048845 

251639.106257 

275441.730057 


229579.796586 

255447.163600 

278925.848153 


815.842923 

849.541882 

875.907645 


2335.246843 

2420.588229 

2470.623727 

m=3,n=2 

3581.204860 

3984.663738 

4371.358784 


225336.052210 

251639.126099 

275441.778797 


229601.353510 

255473.364761 



1186.213960 

1233.926887 

1273.545434 


2714.061057 

2750.657059 

2778.808614 

m=3,n=3 

4709.418836 

5394.586855 

6007.936827 


225336.058337 

251639.160260 

275441.861964 


229637.273129 

255517.021396 

279006.643219 
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Table 4.7: Mean natural frequency for different Web Thickness, lay up sequence (3) 


Mode Shape 

90707i7/0 o /90 o 


Thickness=0.3 

Thickness=0.4 

Thickness=0.5 


iOiji — 1 j . . . , 5 

= 1,... ,5 

= 1, — ,5 

B 

II 

P 

II 

H— 1 

277.546514 

289.136129 

302.529341 

1231.512680 

1289.511493 

1324.384368 

1983.809895 

2183.451137 

2378.925731 

215540.014275 

243925.589135 

269336.015633 

215548.049243 

243935.157969 

269346.913472 

m=l,n=2 

517.546930 

539.802844 

566.371094 

1580.321471 

1594.054759 

1606.586581 

3159.837519 

3634.239119 

4055.834920 

215540.090374 

243925.699134 

269336.174976 

215572.675968 

243964.444321 

269380.283626 

m=l,n=3 

717.441356 

752.611552 

798.086867 

1763.139197 

1780.795611 

1797.331810 

4600.358663 

5308.749981 

5935.792667 

215540.146558 

243925.807910 

269336.361009 

215613.779438 

244013.314688 

269435.962614 

m=2,n=l 

388.709035 

399.130632 

408.742282 

1496.840694 

1627.981165 

1723.276787 

2718.580112 

2844.082161 

2980.744068 

215539.896063 

243925.456174 

269335.870047 

215549.912872 

243937.264623 

269349.240222 
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Table 4.7: (Continued) 


Mode Shape 

90°/0 o /H/0°/90 o 


Thickness=0.3 

Thickness=0.4 

Thickness=0.5 


w„*= 1, — ,5 

cJ x , z = 1 , . . . , 5 

u>„i = l,...,5 


694.780152 

713.121854 

729.172669 


2254.472189 

2296.991193 

2323.287469 

m=2,n=2 

3381.967610 

3826.227849 

4235.475627 


215540.075564 

243925.675222 

269336.145775 


215574.435821 

243966.441477 

269382.493477 


1012.645350 

1039.845355 

1064.352362 


2436.539593 

2453.944665 

2472.913093 

m=2,n=3 

4698.708448 

5403.903875 

6030.752534 


215540.156070 

243925.809434 

269336.358801 


215615.514346 

244015.285662 

269438.144616 


469.595400 

481.156153 

490.541235 


1678.407943 

1828.528562 

1950.694852 

m=3,n=l 

3719.319095 

3825.047251 

3937.480315 


215539.773959 

243925.317069 

269335.717305 


215552.943949 

243940.693154 

269353.028169 


768.248933 

786.281473 

800.237375 


2749.994214 

2913.514825 

3007.015376 

m=3,n=2 

3992.383691 

4315.748387 

4659.862532 


215540.054336 

243925.640868 

269336.104008 


215577.365398 

243969.764511 

269386.169592 


1109.198346 

1135.451955 

1155.587310 


3260.351621 

3303.878696 

3333.763607 

m=3,n=3 

4926.334587 

5604.327792 

6220.720444 


215540.171154 

243925.812158 

269336.355883 


215618.406574 

244018.570372 

269441.780456 
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Figure 4.1: Variance of Natural Frequencies Vs SD of Input Random Variable E 
0 o /90 o /90 o /0°/ir/0 o /90 o /90 o /0 o L/R=1.5, Thickness=0.3 
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Figure 4.10: Variance of Natural Frequencies Vs SD of Input Random Variable G 
0 o /90 o /90 o /0 o /ir/0 o /90 o /90 o /0 o L/R=1.5, Thickness=0.3 
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Figure 4.11: Variance of Natural Frequencies Vs SD of Input Random Variable G 
0° /90 0 /90 0 /0 0 /if/0 0 /90 0 /90°/0 0 L/R=1.0, Thickness=0.3 
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Figure 4.21: Variance of Natural Frequencies Vs SD of Input Random Variable G 
0° /90 o /90 £ ’/0 o /i3 r /0 o /90 o /90 o /0 o L/R=0.5, Thickness=0.3 
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Figure 4.22: Variance of Natural Frequencies Vs SD of Input Random Variable v 12 
0° /90° /90°/0°/iiJr/0 o /90 o /90 o /0 o L/R=1.5, Thickness=0.3 
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Figure 4.27: Variance of Natural Frequencies Vs SD all material properties 
0 0 /90 0 /90 0 /0 0 /i7/0 0 /90 c, /90°/0 0 L/R=0.5, Thickness=0.3 
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Figure 4.28: Natural Frequencies Vs Input Random Variable FEM Result 





Chapter 5 


Conclusion 


The study analysis a honeycomb sandwich laminate composite cylinder. The honey- 
comb core is modified to an equivalent layer of uniform thickness homogeneous material. 
The material properties in the composite laminate as well as the core assumed to be 
random variables. Higher order shear effects have been included in the analysis. The 
laminate is processed for natural frequency by adopting a perturbation approach. The 
mean and various of the natural frequencies for different mode shapes have been eval- 
uated. These are also evaluated by using Monte Carlo simulation with FEM. However 
the sample size generated could not be of representative size because of the large CPU 
time requirement. 

The main conclusion reached on the basis of the above results are: 

1 . In general the first mode natural frequencies have the largest dispersion due to 
variation in the material properties. 

2 . Out of the facing layer material properties £22 has the largest effect of the natural 

frequency variance. 

3. Amongst the material properties for the core E s shows the largest effect on the 
natural frequencies. 

4. The facing layer has a stronger effect on the cylinder natural frequencies than 
the core of the honey comb. 

5. The material property u i2 for the facing layer has the smallest effect on the 
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natural frequencies. The sandwich cylinder shows low sensitivity to G23 and G31 

also. 


5.1 Suggestion for future work 

The present study is limited to linear problem and use only first order perturbation. 
It would be interesting to study nonlinear large deflection problem with higher order 
perturbation procedure. The other related problem in the area are the study of the 
honeycomb sandwich cylinder. The study can continue to different types of cores with 

different thickness. 

The FEM model attempted for the problem resulted in a very large size problem 
with large CPU time requirement. Some rethinking is required to have a smaller 
problem that can be economically handled. 
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A.1 Stiffness Matrix Coefficent 
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-P12P2.P45P54-P12P21P34P43-P12P21P35P53-P12P23P31P43+P12P23P34P41+P12P24P31P43 

»P12P24P33P4X“P12P24P41P55+P12P24P51P45+P12P25P53P31-P12P25P33P51+P12P25P41P54 

-P12P2 5 P51P44-P13P21P32P44-P13P21P32P55+P13P21P34P42-P13P21P35P52+P13P22P31P44 

+ PI3P22P31P55 -P13P22P34P41 “Pl3P22P35pSl + P13P31P44P55 -pi3P3lP5lP45-Pl3P34P4lP55 

•'i-p 1 3P34P45P51+Pl3P35P4lP54-Pl3P35P5lP44-Pl3P24P3lP42+Pl3P24P32P4X-Pl3P25P3lP52 

-4-P13P25P32P51+P14P21P32P43-P14P2XP33P42-P14P21P42P55+P14P21PS2P45-P14P22P31P43 

+P14P22P33P41+P14P22P41P55-P14P22P45P51-P14P31P43P55+P14P45P53P31+PX4P33P41P55 

-PX4P33P45P51-P35P4XP14P53+P14P51P35P43+P14P23P31P42-PX4P23P32P41-P14P25P4XP52 

+PX4P25P42P51+P15P21P32P53+P15P21P33P22+P15P21P42P54-PX5P21P52P44 Pl 5 P22P3lP52 

124 



+P5lP22P33Pl5~Pl5P22P4lP54+Pl5P22P44P51+Pl5P3lP43jP54-Pl5P3lP52P44-Pl5?>33P4lP54 
+Pl5P33P44P51+Pl5P34P4lP52-Pl5P34P43P51+Pl5P23P3lP52-Pl5P23P32P51+Pl5P24P4lP52 
~ P15P24P42P51 - P23P34P42P55 +P23P34P45P52 +P35P42P54P23 
-^5 =PllP22P33P44P55-PllP22P33P45P54-PllP22P34P43P55+PllP22P34P53P45+PllP22P35P43PS4 
-P11P22P35P53P44-P23P11P32P44P55+P45P54P32P23P1I+P23PHP34P42P55-P23PUP34P45P52 

~ P23P UP35P42P54 +P23PllP35P52P44 +P24P1 1P32P43P55 ~ P24Pl 1P32P45P53 — P24PllP33P42P55 

+P24PuP33P45P52+P24PllP35P42P53“P24PllP35P43P52-PllP25P32P43P54+PllP25P32P44P53 
+P33P42P54P1 lP25 -PllP25P33P52P44 -PnP25P34P42P53+PllP25P34P43P52 — P12P21P33P44P55 
+P12P21P33P45P54+P12P21P34P43P55 - Pl2P2lP34P45P53~Pl2P2lP35P43P54+Pl2P2lP35P53P44 
+Pl2P23P3lP43P55 — Pl2P23P3lP45p53~Pl2P23P34P4lP55+Pl2P23P34P5lP45+Pl2P23P35P4lP53 
~ P12P23P35P43P51 - Pl2P24P3lP43P55+Pl2P24P3lP53P45+Pi2P24P33P4lP55— P12P24P33P51P45 
-Pl2P24P35P4lp53+Pl2P24P35P5lP43+Pl2P25P3lP43P54-Pl2P25P3lP44P53 _ Pl2P25P33P4lP54 
~ Pl2P25P33P5lP44+Pl2P25P34P4lP53 — Pl2P25P34P43P51+Pl3P2lP32P44P55~Pl3P2lP32P45P54 
-Pl3P2lP34P42P55+Pl3P2lP34P52P45-Pl3P2lP35P42P54+Pl3P2lP35P52P44-Pl3P22P3lP44P55 
+P13P22P51P45P31+P43P22P34P41P55-P13P22P34P45P51-P13P22P35P41P54+P13P22P35P51P44 
+Pl 3P24P31P42P55 “P 13P24P31P52P45 “ . Pl3P24P32P4lP55+Pl3P24P32P45P51 ~ Pl3P24P35P4lP52 
+P 1 3P24P42P51P35 -P 13P25P31P42P54 +P13P25P3 1P44P52 +P13P25P32P41P54 — P13P25P32P44P51 
-P13P25P34P41.P52+P13P25P42P51P34-P14P21P32P43P55+P14P21P32P45P53+P14P21P33P42P55 
-P14P21P33P52P45-P14P21P35P42P53+P14P21P35P52P43+P14P22P31P43P55-PX4P22P45P53P31 
~P 14P22P33P41P55 +P14P22P33P45P51 +P14P22P35P41P53 -P14P22P51P35P43 -P14P23P31P42P55 
+PhP23P3XP52P45 +P4lP23P32P4lP55 “PX4P23P32P51P45 “P14P23P35P41P52 +P14P23P35P42P51 

+P14P25P31P42P53-P14P25P31P52P43-P14P25P32P41P53+P14P25P32P43P51+P14P25P33P41P52 
-P14P25P33P42P51+P15P21P32P43P54-P15P21P32P53P44-P15P21P33P42P54+P15P21P33P52P44 
+P1SP21P34P42P53-P15P21P34P43P52-P15P22P31P43P54+P15P22P31P52P44+P15P22P33P41P54 
— P15P22P33P44P51 - P34P4lP52Pl5P22+Pl5P22P34P43P51+Pl5P23P3XP42P54 — P15P23P31P52P44 
-P15P23P32P41P54+P15P23P51P44P32+P15P23P34P41P52-P15P23P34P42P51-P15P24P31P42P53 
+P15P24P31P52P43+P15P24P32P41P53-P15P24P32P43P51-P15P24P33P41P52+P15P24P33P42P51 
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